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We review and complete the existing literature on the kinetic theory of spatially homogeneous 
systems with long-range interactions taking collective effects into account. The evolution of the 
system as a whole is described by the Lenard-Balescu equation which is valid in a weak coupling 
approximation. When collective effects are neglected it reduces to the Landau equation and when 
collisions (correlations) are neglected it reduces to the Vlasov equation. The relaxation of a test 
particle in a bath is described by a Fokker-Planck equation involving a diffusion term and a friction 
term. For a thermal bath, the diffusion and friction coefficients are connected by an Einstein relation. 
General expressions of the diffusion and friction coefficients are given, depending on the potential of 
interaction and on the dimension of space. We also discuss the scaling with A'' (number of particles) 
or with A (plasma parameter) of the relaxation time towards statistical equilibrium. Finally, we 
consider the effect of an external stochastic forcing on the evolution of the system and compute the 
corresponding term in the kinetic equation. 



I. INTRODUCTION 



Recently, the dynamics and thermodynamics of systems with long-range interactions has been a subject of consid- 
erable interest in statistical mechanics Some of these systems, perhaps the most important, have been studied 
for a long time. This is the case of Coulombian plasmas, self-gravitating systems and two-dimensional vortices. In 
the last few years, other systems with long-range interactions have been introduced (or revived) and studied. They 
include one and two dimensional plasmas, one and two dimensional self-gravitating systems and toy models such as the 
Hamiltonian mean field (HMF) model and the one dimensional ring model Q ■ Although these models are relatively 
academic, they are of great conceptual interest to better understand the peculiar properties of long-range interacting 
(LRI) systems. Furthermore, they are simpler to simulate numerically and this is certainly valuable to perform very 
accurate comparisons between numerical experiments and theory. Finally, these toy models may find unexpected 
applications in physics. For example, the HMF model presents many features analogous to the Colson-Bonifacio 
Hamiltonian model describing single-pass free electron lasers (FEL). 

The equilibrium statistical mechanics of systems with long-range interactions is now relatively well understood Q ■ 
One important feature of these systems is that they may exhibit ensemble inequivalence and negative specific heats 
first observed in astrophysics (see [B-Q] for reviews) . The out-of-equilibrium properties of these systems are also very 
interesting. In particular, long-range interacting systems may be stuck in non-Boltzmannian quasi stationary states 
(QSS) that persist for a very long time before finally reaching the Boltzmann distribution. These QSSs are stable 
steady states of the Vlasov equation which governs the "collisionless" evolution of the system. In order to understand 
the different timescales involved in the dynamics, it is important to develop a kinetic theory of systems with long-range 
interactions. Since these systems are generically spatially inhomogeneous, the kinetic theory is relatively complicated 
and requires the introduction of angle-action variables [8l-[ll|. However, in some cases (e.g. plasma physics), a stable 
homogeneous phase exists, and, in other cases (e.g. stellar dynamics), it is possible to make a local approximation and 
treat the system as if it were spatially homogeneous. If we restrict ourselves to spatially homogeneous systems, most of 
the work has already been done in the context of plasma physics and stellar dynamics since the kinetic equations that 
govern the "collisional" evolution of a spatially homogeneous plasma or a stellar system, namely the Landau, Lenard- 
Balescu and Fokker-Planck equations, can immediately be transposed to other systems with long-range interactions. 
However, in the details, the diffusion and friction coefficients, and the relaxation time, sensitively depend on the form 
of the potential of interaction and on the dimension of space. It is therefore useful to develop a general formalism 
common to all systems with long-range interactions and, from this general formalism, describe specific models. 

Since the kinetic equations initially introduced in plasma physics and stellar dynamics are now studied by a larger 
community, it may be of interest to regroup in a self-consistent paper the derivation of the basic equations (Vlasov, 
Landau, Lenard-Balescu and Fokker-Planck) to the attention of a wider audience. This is done in this paper (paper I) 
where we extend the derivation given in to a more general class of potentials of interaction in arbitrary dimensions 
of space. For the unity of the subject, we follow their presentation as closely as possible. We also consider the effect 
of an external stochastic forcing on the evolution of the system and compute the corresponding term in the kinetic 
equation. In the second paper (Paper II), we use these general results to determine the explicit expressions of the 
diffusion and friction coefficients, and the relaxation time, for several systems with long-range interactions. Our two 
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papers are therefore essentially pedagogical and aim to synthesize and develop the kinetic theory of systems with 
long-range interactions. They also show what is common to all systems with long-range interactions and what is 
specific to such or such system. 



II. EVOLUTION OF THE SYSTEM AS A WHOLE 



A. Klimontovich and Vlasov equations 

Let us consider a system of material particles in interaction in a d-dimcnsional space. Their dynamics is described 
by the Hamilton equations 

dvi dH dvi dH 

dt ovi dt OTi 

with the Hamiltonian 

^ 1 

■i=l i<j 

where M(|r — r'|) is the potential of interaction and m the individual mass of the particles. For an isolated system, 
the total energy E = H is conserved. We assume that the potential of interaction decays at large distances as 
with 7 < (i. In that case, the potential is said to be long-ranged Long-range potentials include the gravitational 
and the Coulombian potentials, corresponding to 7 = d — 2. 

We introduce the discrete distribution function fd{v,v,t) = ~ ~ '^i(^))- Differentiating this 

expression with respect to time and using the equations of motion (Il|)-([2|). we find that /d(r, v, t) satisfies an equation 
of the form 

dfd , dfd dfd , , 

_+v. — -V<I>.. — =0, (3) 

where 

Mr,t)= J u{\r-r'\)f4r\V,t)dr'dV, (4) 

is the discrete potential produced by /^(r, v,i). This equation is exact and contains the same information as the 
Hamiltonian system ([T|)-(l2]). In plasma physics, it is called the Klimontovich equation [l^ . 

We now introduce a smooth distribution function /(r, v,t) = (/^(r, v,t)) corresponding to an average of /^(r, v,i) 
over a large number of initial conditions. We then write /d = / + <5/ where Sf denotes the fiuctuations around the 
smooth distribution. Similarly, we write = ^ + S^. Substituting these decompositions in Eq. ([3]), we obtain 

df 85 f df 05 f „^ df d5f df d5f „ 

— H + V • — + V • — - - V<i> • — - V$ • — - - V(5<i> • — - V(5$ • — - = 0. (5) 

dt dt dr dr dv dv dv dv 

Taking the average of this equation over the initial conditions, we get 

% + ^-'jf-y^-'i-^-{sm^), (6) 

dt dr ov dv 

where the right hand side can be interpreted as a "collision" term. It docs not really correspond to direct collisions 
between particles but rather to close encounters, or correlations, due to finite N effects. Subtracting this expression 
from Eq. ([5]), we obtain an equation for the perturbation 

^ + V . ^ - V<i> . - . ^ = ^ . iSfW5<,) - A . ^SfW5^) . (7) 
dt dr ov ov ov ov 

These equations are still exact since no approximation has been made so far. 

Using dimensional analysis in Eq. ([2]), we find that the kinetic energy and the potential energy are comparable 
if Nmv^ ~ N'^m?u. Therefore, the energy, the temperature and the dynamical time scale like E ~ N'^m'^u, (3~^ ~ 
mv^ ~ Nm^u and to ~ R/v ^ R/^/Nmu (where R ^ V^l'^ is the system size). We now consider the thermodynamic 
limit N — > -|-oo in such a way that the normalized energy e = E l{N'^vr?u\ the normalized temperature t] = jSNrn^u 



3 



and the normalized dynamical time td ~ toV Nmu/ R are of order unity. By a suitable normalization of the physical 
quantities, it corresponds to — +00 with V ^ 1, m ^ 1, E ^ N, (5 ^ 1, to ^ and u ~ This is equivalent 

to putting in front of the potential energy in the Hamiltonian (Kac prescription) [l^l- Since the strength of 
the potential goes to zero as — >■ +00, we are dealing with weak long-range potentials. In this scaling, the energy 
is extensive and the temperature intensive. Furthermore, f/N ~ 1, $ ^ 1, 6f/N ~ 1/VN and (5$ ~ 1/\/N. We 
can therefore consider an expansion of the equations of the problem in terms of the small parameter 1 /N which is a 
measure of the "graininess" of the distribution (as approaches zero, coUisional effects disappear) [g^l- Considering 
Eq. ([5]), we sec that the collision term (r.h.s.) scales like 1/N as compared to the advective term (l.h.s.). Therefore, 
for N +00, Eq. ^ reduces to the Vlasov equation 

^+v-|^-V$.|^=0, <i>(r, t)^ j u{\v - r'|)/(r', v', i) dv'dV. (8) 

The Vlasov equation is obtained when "collisions" (more properly granular effects, discreteness effects or finite N 
corrections) between particles are neglected. The Vlasov equation can also be obtained from the Liouville equation, 
writing the BBGKY hierarchy and neglecting correlations between particles. This is valid, for a fixed interval of time, 
when N — > +00. In that case, the mean field approximation becomes exact and the TV-body distribution function can 
be factorized in a product of N one-body distributions, resulting in Eq. (|5]). The "Vlasov equation", which is also 
called the colhsionless Boltzmann equation (CBE), was introduced by Jeans in stellar dynamics and by Vlasov 
in plasma physics. In fact, it is valid for a large class of systems with long-range interactions The exactness 
of the Vlasov equation when — > +C)0 was proven rigorously by Braun & Hepp [l6[ for non singular potentials. 

Starting from an unsteady or unstable initial condition, the Vlasov equation, coupled to a long-range potential 
of interaction, is known to develop a complicated mixing process leading to the formation of a quasistationary state 
(QSS) on a coarse-grained scale. This process, which is purely colhsionless and driven by mean field effects, is known as 
violent relaxation. It takes place on a very short timescale, of the order of a few dynamical times ~ 1. A statistical 
theory of the Vlasov equation has been developed by Lynden-Bell in astrophysics to predict the QSS that results 
from violent relaxation, assuming ergodicity (or, at least, efficient mixing). A similar statistical theory has been 
developed by Miller [l^ and Robert & Sommeria [l^ for the Euler equation in two-dimensional (2D) hydrodynamics. 
The process of violent relaxation towards a non-Boltzmannian QSS has also been extensively studied for the HMF 
model Q. In some cases, violent relaxation is incomplete due to lack of ergodicity or efficient mixing [20| . 

On longer timescales, "collisions" (more properly correlations) between particles develop and the system deviates 
from the Vlasov dynamics. For t +00, we expect that the system will reach a statistical equilibrium state. In 
the thermodynamic limit N — > +00, the mean field approximation is exact [2l|. The statistical equilibrium state 
in the microcanonical ensemble is obtained by maximizing the Boltzmann entropy S = — J {f /m)hi{f /m) drdv at 
fixed mass M = j f drdv and energy E = (1/2) / fv^ drdv + (1/2) / p$ dr where p{r, t) ~ J /(r, v, t) dv is the mass 
density. This variational principle determines the most probable macrostate /(r, v). Writing the first order variations 
as 6S — /3SE — a5M = where 13 and a are Lagrange multipliers, we obtain the mean field Boltzmann distribution 

/(r,v) = Ae-^™('^+*('-», <i>{r) ^ J u{\r-r'\)p{r')dr'. (9) 

For the Newtonian interaction (gravity) , the equilibrium potential is solution of the Boltzmann- Poisson equation. In 
the case of plasmas with a neutralizing background, the statistical equilibrium state is spatially homogeneous and 
Eq. (O reduces to the Maxwell distribution. We must remember, however, that the statistical theory is based on an 
assumption of ergodicity (or efficient mixing) and on the postulate that all the accessible microstates are equiprobable. 
There is no guarantee that this is always true for systems with long-range interactions. In order to prove that the 
system reaches the Boltzmann distribution (|9]) , and in order to determine the relaxation time (in particular its scaling 
with the number N of particles), we must develop a kinetic theory of systems with long-range interactions 0, |^, [22j . 

B. Landau and Lenard-Balescu equations 

As we have previously indicated, the "collisions" between particles can be neglected for times shorter than Ntu 
where i/j is the dynamical time. If we want to describe the evolution of the system on a longer timescale, we must 
take finite N corrections into account. At the level 1/A^, we can neglect the quadratic term on the right hand side of 
Eq. ([7]). Indeed, the left hand side is of order 1/Vn and the right hand side is of order 1/A^ <^ 1/\/N . We therefore 
obtain the set of coupled equations 



(10) 
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dSf dSf ^ df dSf , , 

at or c/v ov 

They form the starting point of the quasihnear theory whieh is vahd in a weak couphng approximation (23j . It can be 
shown that these equations describe the evolution of the system under the effect of two-body colhsions (higher order 
correlations are neglected). If we restrict ourselves to spatially homogeneous distributions, the field — V$ vanishes. 
In that case, Eqs. (fTO|) - (fTT|) reduce to the coupled equations 



These equations are valid as long as the spatially homogeneous distribution function is Vlasov stable so that it evolves 
under the sole effect of collisions. We shall assume that the fluctuations evolve rapidly compared to the transport 
time scale, so that time variation of / and $ can be neglected in the calculation of the collision term (Bogoliubov 
ansatz). Therefore, for the purpose of solving Eq. (|13p and obtaining the correlation function (SfS^), we shall regard 
/(v) as constant in time. Actually, the time dependence of / is given by Eq. p^ . Only the fluctuation term ((5/(5$) 
drives / and then its evolution is much slower than the evolution of 6f. Indeed, / changes on a timescale of the order 
Nto (or larger) while ((5/(5$) relaxes, by Landau damping, to its asymptotic form on a much shorter timescale, of 
the order of the dynamical time t]j. After the correlation function has been obtained as a functional of /, the time 
dependence of / may be reinserted in the kinetic equation. This is an adiabatic hypothesis which is valid for iV ^ 1. 
With this approximation, Eqs. p^ - ([T5|) can be solved with the aid of Fourier-Laplace transforms and the collision 
term can be explicitly calculated. This is the approach presented in psj to derive the Lcnard-Balcscu equation. 
The Fourier-Laplace transform of the fluctuations of the distribution function 6f is defined by 



dr 



{2n) 



+ 00 



(5/(k,v,w)= / ^— ^ / dte-*(''-"*)(5/(r,v,t). (14) 







This expression for the Laplace transform is valid for Im(w) sufficiently large. For the remaining part of the complex 
Lj plane, it is defined by an analytic continuation. The inverse transform is 

5/(r,v,t) = J (ik^^e'('' '-"*)5/(k,v,c.), (15) 

where the Laplace contour C in the complex lo plane must pass above all poles of the integrand. Similar expressions 
hold for the fluctuations of the potential i5$. We note that, for periodic potentials, the integral over k is replaced by 
a discrete summation over the different modes. If we take the Fourier-Laplace transform of Eq. (|T3)) . we find that 

df ~ 

- (5/(k, V, 0) - iuj (5/(k, V, w) -I- ik • V (5/(k, v, w) - ik ~ (5$(k, lu) ^ 0, (16) 
where the first term is the spatial Fourier transform of the initial value 

/dr 
^ e-^'^-Sf {r,v,0). (17) 

The foregoing equation can be rewritten 

(5/(k,v,..) = JlJLj|.(k,..) + (18) 
k • V — i(k • V — w) 

where the first term on the right hand side corresponds to "collective effects" and the second term is related to the 
initial condition. The fluctuations of the potential are related to the fluctuations of the distribution function by a 
convolution 

(5$(r,t)= J u{\r-r'\)Sfir',v',t)dr'dV. (19) 
Taking the Fourier-Laplace transform of this equation, we obtain 

(5$(k,w) = (27r)''u(fc) J (5/(k,v,w)dv. (20) 
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Substituting Eq. (|T8|) in Eq. (l20l) . we find that the Fourier-Laplace transform of the fluctuations of the potential is 
related to the initial condition by 

SHKu.)^i2nr^fd.^I^^ (21) 
e(k,w)J i(k-v-a;) 

where the dielectric function is defined by 

r k- ^ 

e(k,a;) = 1 - (27r)''u(fc) / ^ dv. (22) 

/ k • V — cj 



The Fourier-Laplace transform of the fluctuations of the distribution function is then given by Eq. with Eq. (|2ip . 
In these equations, and in the following, / stands for /(v) and /' stands for /(v'). The dispersion relation associated 
with the linearized Vlasov equation corresponds to e(k, w) = (see Appendix [B|) . If the system is Vlasov stable, 
e(k, oj) does not vanish for any real uj, so Eq. (|2ip is well-defined. If collective effects were neglected in Eq. p^ . we 
would obtain Eq. (PT|) with e(k, w) = 1. This shows that, due to collective effects, the bare potential of interaction 
u{k) is replaced by a "dressed" potential u(A:)/e(k, cj) taking into account the polarization of the medium. 

We can use these expressions to compute the collision term appearing on the right hand side of Eq. (|12p . One has 



((5/V,5$) = / dk / — / dk' / — zk'e'('' '-"*)e^('^' "-"'*)((5/(k,v,L^)<5l>(k',w')>- 
J Jc J Jc 27r 

Using Eq. we find that 



(23) 



(W,v,.)5^(k',.')) - ^^{5^K.)SH^',.')) + (W,y,O)^^(k>0) ^ (24) 

k • V — i(k • V — w) 

The first term corresponds to the self-correlation of the potential, while the second term corresponds to the correlations 
between the fluctuations of the potential and the fluctuations of the distribution function at time t = 0. Let us consider 
these two terms separately. 
From Eq. (|2ip . we obtain 

u^^v' '^\ fo ^2d^(fcMfc:^ f, , , (J/(k,v,0)^/(k^v^O)) 
((5$(k,w)5$(k ,w )) = -{2tt) / dvdv — — — — — . (25) 



e(k, uj)e(k', lu') J (k • v - w)(k' • v' - w') 

Using the expression of the auto-correlation of the fluctuations at t = given by (see Appendix |^ : 

(<5/(k, V, 0)<5/(k', v', 0)) = j^^Sik + k')<5(v - v')m/(v), (26) 

we find that 

{S^k,u;)S^{k',u')) = {2nrm ,J {k + k') f dv ^[j^ (27) 

e(k, w)e(— k, w') J (k • v — c<;)(k • v + cj') 

Considering only the contributions that do not decay in time, it can be shown [2^ that [(k • v — w)(k • v + w')]^^ 
can be substituted by (27r)^(5(w -I- ijj')S{k ■ v — w). Then, using the property e(— k, —u) = e(k, w)*, one finds that the 
correlations of the fiuctuations of the potential are given by 

(<5$(k, w)(5$(k', w')> - (27r)'^+^m , A(k + k')(5(w + uj') / 5(k • v - w)/(v) dv. (28) 

|e(k,w)P 7 

Similarly, one finds that the second term on the right hand side of Eq. (|24p is given by 

{Sf{k,^,0)S^ik',Lo')) 2 ^(^ArV^VU^ ^ 'U^V \ (OQ^ 

^ = 27r TO— — -5 k + k )S{uj + UJ (5 k • v - a;)/(v . (29 

i(k ■ V — UJ) e(k,uj') 

From Eq. (|28p , we get the contribution to Eq. of the first term of Eq. (IMl) . It is given by 



(5/V<5$){ = -»(27r)'^+iTO fdk [ rfVfc, /'^^ , 

^ ^ ^ J Jc2ttJ *k-v-w|e(k,L^) 



/(v')5(k.v'-u;). (30) 
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Using the Landau prescription w — > w + iO^ and the Plemclj formula, 

^-r[^T^^6i.), (31) 

where V denotes the principal value, we can replace l/(k • v — — iO+) by +i7r(5(k • v — lu). Then, integrating on w, 
we obtain 

{SfVS<^>)l = n{2nrm J dkdV k,k,j-^^^^^^6[k ■ - v')]f{v')^^{v). (32) 

From Eq. (|29p . we get the contribution to Eq. ((23| of the second term of Eq. ((24|) . It is given by 

{5fW6^)'J = m I dkk,j-^^^^^Ime{k,k-^)f{v). (33) 

Using the Landau prescription cj — > cj + iO^ and the Plemelj formula pip , the imaginary part of the dielectric function 
(1221) is 

Im e(k, w) = -tt(2tt fu(k) [ k ■ ^5{k • v - w) dv. (34) 
Substituting this expression in Eq. ()33p . we obtain 

Finally, regrouping Eqs. ()12p . (|32|) and (|35|) . we end up on the kinetic equation 

This equation, taking collective effects into account, is the Lenard-Balescu equation of plasma physics. It was originally 
derived in a very different manner [2^. [25j and, in fact, several derivations exist in the literature. 
If we neglect collective effects and take |e(k, k • v)| = 1, Eq. ([M]) becomes 

^ = ^(2^)'m A J dk dV hk,u{kf5[k ■ (v - v')] - ^) /(v, 0/(v', t). (37) 

The integral over k can be done explicitly (see, e.g., [^) and we obtain 

df d f , .w'^Sii — WiWi f d d \ . , ^ 

it = ^^^d^ j - a^j ^(-'^)^(-' 

where w = v — v' is the relative velocity and Kd is a constant with value = Sn^m J^°° k^u{k)'^ dk in d = 3 and 

K2 = Stt^to /q^°° k^u{ky dk in d = 2. This kinetic equation, that neglects collective effects, is the Landau equation 
[26| . We note that the structure of this kinetic equation does not depend on the specific form of potential of interaction. 
The potential of interaction only appears in factor of the integral in a constant Kd that determines the relaxation 
time. For a 3D plasma, using (2Tr)^u{k) = Ane'^ /m?k^ , we get = (27re''/m'^) In A where In A = J^°° dk/k is a 
Coulomb logarithm that has to be regularized with appropriate cut-offs. The large-scale cut-off is the Debye length 
A|) ^ ksT/ne^ and the small-scale cut-off is the Landau length Al ^ e^/mvf^ ~ e^/ksT ~ 1/nXjj (see Paper II 
for more details). Therefore A ^ Xd/^l ^ ^^A'p 3> 1 represents the number of electrons in the Debye sphere. For 
a 2D plasma, using (27r)^M(fc) = 27re^/m^fc^ and introducing a large-scale cut-off at the Debye length, we obtain 
A'2 = /m^kD- The validity of these heuristic procedures will be further discussed in Paper II. Returning to Eqs. 
([55)1 and ((571) . note that collective effects can be taken into account simply by replacing the bare potential u{k) 
in the Landau equation by a "dressed" potential Udressed{k) = u(fc)/|e(k, k • v)|, including the dielectric function, 
without changing the overall structure of the kinetic equation. Physically, this means that the particles are "dressed" 
by their polarization cloud. In plasma physics, collective effects are important because they account for screening 
effects and regularize, at the scale of the Debye length, the logarithmic divergence that occurs in the Landau equation. 
This avoids the introduction of ad hoc cut-offs. 
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Systems with long-range interactions are gene rically spatially inhomogeneous. If wc neglect collective effects, we 
can derive a generalized Landau equation [ol l27l [28j : 



dl 
dt 



df df d 

W ■ — V$ • — 

dv N dv 



dr J drldvlF^'{l ^ 0)g{t,t - t) 



/(r,v,t-T) — (ri,vi,t~T), 
m 



(39) 



that is valid for systems that arc not necessarily spatially homogeneous, and not necessarily Markovian [6J]. For 
spatially homogeneous Markovian systems, we recover the Landau equation ([37]) after simple calculations Equation 
(p9| shows that the Landau equation has the structure of a Fokker-Planck equation and that the coefficients of diffusion 
and drift are given by generalized Kubo formulae. If we assume that the system is in a quasi stationary state (QSS) of 
the Vlasov equation slowly evolving under the effect of "collisions" , we can derive a spatially inhomogeneous Landau 
or Lenard-Balescu equation written in angle-action variables |8l-[ll|. In certain cases, e.g. for 3D self-gravitating 
systems, we can make a local approximation. In that case, the kinetic equation (j36p is replaced by 



dt 



dl 
dr 



dl 
dv 



_d_ 

dvi 



dkdv' kikj 



u{kf 



|e(k,k. 



,<5[k.(v-v')] 



_a d_ 

dvi dv'. 



/(r,v,i)/(r,v',t). (40) 



In this equation, the effects of spatial inhomogeneity are kept only in the advection (Vlasov) term, while the collision 
term is calculated as if the system were spatially homogeneous. If we neglect collective effects, taking |e(k, k • v) | = I, 
this equation can be rewritten 



df df „^ df d I , ,w 

dt dr dv dvi 



TT-'TTt) /(r,v,t)/(r,v',i). 

dvj dv'j J 



(41) 



For a 3D self-gravitating system, using (27r)'^u(fc) = —AnG/k^, wc get = 2711710^111 A where InA = J^°° dk/k is 
the Coulomb factor that has to be regularized with appropriate cut-offs. The large-scale cut-off is the Jeans length 
ksT / Grn^n (which is of the order of the system size R) and the small-scale cut-off is the gravitational Landau 



>2 



length Al ~ Gm/i 



Gm^/ksT - 1/nAj. Therefore A - Aj/Al 



^ 1 represents the number of 



stars in the cluster. The Vlasov-Landau kinetic equation (|4ip is the fundamental kinetic equation of stellar dynamics 
For a 2D stellar system, using (27r)^'u(fc) = ~2TTG/k^ and introducing a cut-off at the Jeans length, we obtain 



m 

K2 



2'kG m/kj, but the local approximation is not justified in that case (see Paper II). 



C. The relaxation time of the system as a whole 

The Lenard-Balescu equation ([55)1 is valid at the order 1/A^ so it describes the "collisional" evolution of the system on 
a timescale ~ Nto. This kinetic equation conserves the mass M and the energy E (which reduces to the kinetic energy 
for a spatially homogeneous system) . It also monotonically increases the Boltzmann entropy S (iJ-theorem) [l^ . The 
collisional evolution of the system is due to a condition of resonance. This condition of resonance, encapsulated in 
the 5-function, corresponds to k • v' = k • v with v' ^ v. 

When d > 1, the condition of resonance can always be satisfied and the Boltzmann distribution ^ is the unique 
steady state of the Lenard-Balescu equation. Therefore, under the effect of "collisions", the system reaches the 
Boltzmann distribution on a relaxation time 

t/j ^ NtD, {d > 1 homogeneous). (42) 

For 3D plasmas, we get from Eq. ([55)1 an estimate of the relaxation time as tji ^ m'^v'^/[ne^\n{mv^/e^n^^^)] ^ 
A^/^A^^m^/^/eln A, leading to tfj ^ (A/lnA)t£). For 2D plasmas, we obtain tfj ^ im?v^^kDlne^ ^ m^/^v^/n-^/^e"^ 
A^/^A/jm^/^/e, leading to tj^ ^ Atjj. In both cases, A = nA^ represents the number of electrons in a Dcbye sphere 
and tjj ^ ujp^, where wp ^ (ne? jmYl"^ ^ Vm/\D is the plasma pulsation, is a measure of the dynamical time. For 
3D stellar systems, we get from Eq. (PT|) an estimate of the relaxation time as tp ^ v'^/[nm'^G^ ln(z;^/Gmn^/^)] ~ 
^1/2^^3/2^^1/2^1/2 leading to tp {N/lnN)tD. For 2D stellar systems, we obtain tR ~ vl^kj /nw?G'^ ~ 
w^/n-^/^m'^/^G^/^ '-^ A^^/^A,//m^/^G^/^ ~ Nt^- In both cases, A = nAj ~ N represents the number of stars in the 
cluster and tjj ^ '-^g^^ where ujq ^ (Gp)^^^ ^ v,„/Aj is the gravitational pulsation, is a measure of the dynamical 
time. We have just given the scaling of the relaxation time as the prefactor depends on the precise way the relaxation 
time is defined. 
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For spatially homogeneous one dimensional systems, there is no resonance and the Lenard-Balescu collision term 
vanishes. Indeed, the kinetic equation p6p reduces to 

so the distribution function does not evolve at all on a timescale ^ Nto. Therefore, the kinetic theory predicts no 
thermalization to a Maxwellian at first order in The maxwcllization is at least a second order effect in and 

consequently a very slow process. The relaxation time satisfies 

tfi > Nto (d = 1 homogeneous). (44) 

This result has been known for a long time in plasma physics [3^ . [ssj and was rediscovered recently in the context 
of the HMF model [36[ [3^ . Since the relaxation process is due to more complex correlations than simply two-body 
collisions, we have to develop the kinetic theory at higher orders (taking three-body, four-body,... correlation functions 
into account) in order to obtain the relaxation time. If the collision term does not vanish at the next order of the 
expansion in powers of 1/iV, the kinetic theory would imply a relaxation time of the order N^tu. This quadratic 
scaling, conjectured in is observed for ID plasmas that are spatially homogeneous [H, In that case, the 
relaxation is caused by three-body correlations. However, the problem could be more complicated and yield a larger 
relaxation time like e^to, as observed in [i^l for the permanently homogeneous HMF model (but other authors |4]| 
find the more natural A^^ scaling). In fact, it is not even granted that the system will ever relax towards statistical 
equilibrium; the evolution may be non-crgodic (the mixing by "collisions" may be inefficient). 

For spatially inhomogeneous systems, since there are potentially more resonances (as can be seen from the Lenard- 
Balescu equation written in angle-action variables |8l4llj). the relaxation time could be reduced and approach the 
natural scaling 

tji ^ Nto (inhomogeneous), (45) 
predicted by the first order kinetic theory [1, Q. This linear scaling is observed for ID self-gravitating systems that 



arc spatially inhomogeneous |42h45| and for the permanently inhomogeneous HMF model [46|. However, very little 
is known concerning the properties of the spatially inhomogeneous Lenard-Balescu equation and its convergence (or 
not) towards the Boltzmann distribution. It could approach the Boltzmann distribution (since entropy increases) 
without reaching it exactly if, at some point, the condition of resonance cannot be satisfied anymore. 

Finally, if we consider a ID system that is initially spatially homogeneous but, slowly evolving under the effect of 
collisions, becomes Vlasov unstable and undergoes a dynamical phase transition [47| making it spatially inhomoge- 
neous, wc expect a relaxation time scaling like with 1 < i5 < 2, intermediate between the scalings pij) and of 
permanently homogeneous and permanently inhomogeneous systems. Such a scaling, with 8 = 1.7, has been observed 
in [4^ for the HMF model experiencing a dynamical phase transition. 

The same type of results are found for 2D point vortices [1^ . 

III. STOCHASTIC PROCESS OF A TEST PARTICLE: DIFFUSION AND FRICTION 

A. The Fokker-Planck equation 

In the previous section, we have studied the evolution of the system as a whole. We now consider the relaxation of a 
test particle in a bath of field particles with a steady distribution function /(v). We assume that the field particles are 
either (i) at statistical equilibrium with the Boltzmann distribution (thermal bath), in which case their distribution 
does not change at all, or (ii) in the case of one-dimensional systems, in any stable steady state of the Vlasov equation 
(as we have just seen, in c? = 1, this profile does not change on a timescale of order Ntjj). We assume that the 
test particle has an initial velocity Vq and we study how it progressively acquires the distribution of the bath due to 
"collisions" with the field particles. As we shall see, the test particle has a stochastic motion and the evolution of 
the distribution function P(v, t), the probability density that the test particle has a velocity v at time t, is governed 
by a Fokker-Planck equation involving a diffusion term and a friction term that can be analytically obtained. The 
Fokker-Planck equation may then be solved with the initial condition P(v, to) = S{v — vq) to yield P(v, t). In plasma 
physics, this test particle approach has been developed by Hubbard [50| who properly took collective effects into 
account. The derivation that we give below is slightly different and follows [3] closely. 

The equations of motion of the test particle are 

|=v, 'i--ySHr,t). (46) 
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When TV — )• +00, the fluctuations of the potential can be neglected and the test particle follows a linear trajectory 
at constant velocity v. At the order the fluctuations (5$ of the potential must be taken into account and the 

particle has a stochastic motion. Equations (|46| can be formally integrated into 

r(t) = r+ / wit')dt', (47) 
Jo 

v(t)=v-/" W6<i>ir{t'),t')dt', (48) 

where we have assumed that, initially, the test particle is at (r,v). Note that the "initial" time considered here does 
not necessarily coincide with the original time to mentioned above. Since the fluctuation (5$ of the potential is a small 
quantity, the foregoing equations can be solved perturbatively. At the order l/iV, which corresponds to quadratic 
order in (5$, we get for the velocity 

n,{t)^u.- j^dt'^£-{v + .t'X)+ j^dt' dt" dt"'£^^{r + vt\t')^{r + .r,t"'). (49) 

As the changes in the velocity are small, the dynamics of the test particle can be represented by a stochastic process 
governed by a Fokker-Planck equation [51[. If we denote by P{v,t) the density probability that the test particle has 
a velocity v at time the general form of this equation is 

dP 52 d 

The diffusion tensor and the friction force are defined by 

Aj(v) = lim ^(Mt) - v,)iv,it) - V,)), (51) 

t— S-+00 Zt 

A,;(v)= lim l{v^{t)-v^). (52) 

In writing these limits, we have implicitly assumed that the time t is long with respect to the fluctuation time but 
short with respect to the relaxation time (of order Nto), so the expression p9)) can be used to evaluate Eqs. (|5T|) 
and (j52p . As shown in our previous papers |32| . [52j . it is relevant to rewrite the Fokker-Planck equation in the 
alternative form 



dP d f„ dP 

The total friction is 



^-^\D,~-PFr']. (53) 



A,=Fr' + ^, (54) 

where Fp^i is the friction due to the polarization j6^ . while the second term is due to the variations of the diffusion 
tensor with the velocity v. As we shall see, this decomposition arises naturally in the following analysis. The two 
expressions ([50]) and ([55]) have their own interest. The expression ([501) where the diffusion tensor is placed after the 
second derivative {DP) involves the total (dynamical) friction A and the expression (|53)) where the diffusion tensor 
is placed between the derivatives dDdP isolates the friction by polarization Fpoi. We shall see in Sec. IIII El that this 
second form is directly related to the Lenard-Balescu equation (|36p . It has therefore a clear physical interpretation. 
We shall now calculate the diffusion tensor and the friction term from Eqs. ([51]) and ([5^ . using the results of Sec. 
Ill Bl that allow us to take collective effects into account. 



B. The diffusion tensor 



We first compute the diffusion tensor defined by Eq. ([5T|) . Using Eq. ([^0)) . we see that it is given, at the order 
1/N, by 

1 /■* /■* / 85^ dS^ 
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By the inverse Fourier-Laplace transform, we have 

/f5(r + vt',0f^(r + vt",^")W- /dk fdk' [ ^fc,fc'e»>^(r+vO,-..*'e^k',r+v*")^-..'*" 

\ ^x^ dxj I J Jc 27r 7 Jc 27r 

x((5l>(k,a;)5l>(k',w')). (56) 
Substituting Eq. (|28p in Eq. (|56p . and carrying out the integrals over k', w' and oj, we end up with the result 

This expression shows that the correlation function appearing in Eq. (|55p is an even function of t' — t" . Using the 
identity 

dt' f dt" fit' -t") = 2 f dt' f dt" f{t' -t")^2 [ ds{t-s)f{s), (58) 
"'0 "'0 Jo Jo 

we find, for t — > +oo, that 

A.^/ ( — (r,0) — (r + v.,.))d. (59) 

This is the Kubo formula for our problem. Replacing the correlation function by its expression (|57p . we get 

A, - {27rfm ds f dkdv' fc,fc,e*'^ (^-^>— 4^^L/(v'). (60) 
Jo J |e(k,k-v')r 

Making the change of variables s — > — s and k — > — k, we see that we can replace Jf^°° ds by (1/2) ds in Eq. (|60p . 
Then, using the identity 

/'+°° dt 

.M = /^e.".i, (61) 

we obtain the final expression 

A, = 7r(27r)'^m J dkdV k,k,j-^^^^^S[k ■ (v - v')]/(v'). (62) 

C. The friction due to the polarization 

We now compute the friction force defined by Eq. (|52p . We need to keep terms up to order 1/A^. From Eq. (|49p . 
the first term to compute is 

1 /■* 

A/ = — / dt' (V(5$(r + vt', t')) . (63) 
i Jo 

By the inverse Fourier-Laplace transform, we have 

(V(5$(r + vi',t')) = t f dk f — ke''"(''+^*')e-'"*'((5$(k,cj)). (64) 
Using Eq. ([201), we find that 



(,4(k,.,).,2,,^j|^L/.v'<^«!^. ,65, 

e(k,a;)7 i(k-v'-w) 

Now, using the fact that the test particle is initially located in (r, v), so that {6f{r', v', 0)) ~ 7Jid{r' — r)(5(v' — v), we 
obtain from Eq. p7)) the result 

(5/(k,v',0)) = ^me-"'-5(v' - v). (66) 
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Substituting these expressions in Eq. (|64)) . we get 



(V(5$(r + vt',0> / dk I ^ ke*(k v-^)t'^(^^ 1 . (57) 

' ' 2tt e(k, w) k • V — cj 



Therefore, the friction term (j63p is given by 

A, = -ml f dt' [dkfp ke'(— ^ . (68) 
tj,, J Jc27T e{k,uj)k-v-oj ^ ^ 

We now use the Landau prescription lj ^ lj + iO+ and the Plemelj formula ([31]) to evaluate the integral over to. The 
term corresponding to the imaginary part in the Plemelj formula is 



A^"^ = -imn- / dt' / dk / _ keK''-^-'^^* — ^<5(k • v - w). (69) 
t Jo J J '2'^ e(k,w) 



Integrating over lo and i', we obtain 



Al")=-^/dkk^^^i^ImKk,k.v)]. (70) 



The term corresponding to the real part in the Plemelj formula is 



Af) = -„^i tdt' [dkP ^""^ke-C— ' . (71) 

Integrating over t' , we can convert this expression to the form 

Af^ ^ im f dkV r — k i ^i{isin[(k-v-a;)i]+cos[(k-v-w)t]-l}. (72) 

7 J-oo 27r e(k,w) (k- v-w)^ t 



For t — > +00, using the identity 

1 — cos(te) , 
.1™ 7^^=^-5(2^), (73) 

and integrating over w, we find that A^**' is given by Eq. ([70| . just like A^"'. Therefore, writing A/ = A'f^ + Af' — 
2A'f\ we obtain 

Ai = -m fdkk "f^^ Im [e(k, k • v)] . (74) 
7 |e(k,k-v)|^ 

Finally, using Eq. (p4|. we find that 

A[ = 7r(27r)'^m J dkdV k,k,j-^^^^^6[k ■ (v - ^')]^{^')- (75) 

As wc shall see. A/ corresponds to the friction due to the polarization denoted Fpoi in Eq. ((M)) . 

Remark: We can obtain Eq. (j75p in a more direct manner from Eq. (|67p by using the contour of integration shown 
in Fig. 9 of [1^. In that case, the integral over w is just — 27ri times the sum of the residues at the poles of the 
integrand in Eq. (|S7)) . The poles corresponding to the zeros of the dielectric function e(k, w) give a contribution that 
rapidly decays with time since Im(a;) < (the system is Vlasov stable). Keeping only the contribution of the pole 
a; = k • V that does not decay in time we obtain 

(V5$(r + vi',0) [ dkk— ^4^- (76) 

J e(k,k-v) 

Substituting this result in Eq. we get Eq. then Eq. ([75]) . 
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D. The part of the friction due to the inhomogeneity of the diffusion coefficient 

In the evaluation of the total friction, at the order the second term to compute is 

By the inverse Fourier-Laplace transform, we have 

xe''*''('"+^*"')e-*'^'*"'((5l>(k,a;)5l>(k',w')). (78) 



Substituting Eq. (|28p in Eq. (|78p . and carrying out the integrals over k', w' and w, we end up with the result 

This expression shows that the correlation function appearing in Eq. (|77l) is an odd function of t' — t'" . Using the 
identity 



we find that 



/* dt" f dt'" fit' - t'") = f dt'" it' - t"')f{t' - t'"), 
Jo Jo Jo 



A{^ = ^{2nYm- f dt' f dt'" [ dkdV hk\t' ~ t"')e^<^--')it' -t'") J^W' 



(80) 



t Jo Jo J ^ ' |e(k,k.v')| 

This can be rewritten 

- '^"'"'11- i' / wi^^F^'^''' 

Since the integrand only depends on t' — t'" , using the identity we obtain for t +oo. 

Therefore, similarly to the Kubo formula ([55]) for the diffusion coefhcient Dij, the friction Ajj can be written 

Al^= (|-(r,0)£^(r + v.,.)).d.. (84) 

Making the change of variables s — > — s and k — k, we see that we can replace ds by (1/2) ds in Eq. (|83p . 
Then, using the identity (|6ip. we obtain the expression 

A" = A2nr^n^^ J dkdV k,k,j-^^^^^S[k ■ (v - v')]/(v'). (85) 
This expression involves the diffusion tensor (p^ . so that finally 

A" = m 

Remark: we can obtain this identity in a more direct manner by taking the partial derivative of Eq. (j59p with 
respect to Vj (using the fact that Dij = Dji) and comparing with Eq. 
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E. Connection between the Lenard-Balescu equation (|36|) and the Folcker-Planclc equation (|53|) 

We have established that the diffusion tensor and the friction force are given by 

A, = ni2nrm J dkdV k,k,j-^^^^S[k ■ (v - v')]/(v'). (87) 

A. = .i2.rm I dkdV hk.j^^^^Sik . (V - + f^. (88) 

Comparing Eq. (|88p with Eq. ((5^ . we see that the friction by polarization is 

Fr'=n{2nrm J dkdV hk.j-^^^^^Sik ■ - (89) 

On the other hand, using an integration by parts in the first term of Eq. (j88p . the total friction can be written 

A. . ni2nrm I dkdV k.k,fiV) (A _ Sik . (V - v')]^^^§L. (90) 

The foregoing expressions of the diffusion tensor and friction force, taking collective effects into account, were first 
derived by Hubbard [soj . 

Using Eqs. (|87p and we find that the Fokkcr-Planck equation ([55)1 becomes 

If the system is spatially inhomogeneous, but if a local approximation can be justified, the Fokker-Planck equation 
(|9T|) is replaced by 

|:+.,|e_V*,.).|^ ^ ,(2„',„|- /...v'...,^^ifllp*lk.(v-V)l (A - _!_) P(.v,,)/(.v'). (92, 

where $(r) is the static potential produced by /(r, v). The form of Eq. ([?T|) is very similar to the form of Eq. (|36p . 
This shows that the Fokker-Planck equation (|91l) . with the diffusion coefficient ([57)1 and the friction ([SH]) . can be 
directly obtained from the Lenard-Balescu equation psp by replacing the time dependent distribution /(v',t) by the 
static distribution /(v') of the bath. This procedure transforms an integro-differential equation p6p into a differential 
equation (|9ip . Although natural, the rigorous justification of this procedure requires the detailed calculation that we 
have given here. In fact, we can understand this result in the following manner. Equations (j36p and (j9ip govern the 
evolution of the distribution function of a test particle (described by the coordinate v) interacting with field particles 
(described by the running coordinate v'). In Eq. ([M]) . all the particles are equivalent so the distribution of the field 
particles f{v',t) changes with time exactly like the distribution of the test particle f{v,t). In Eq. the test 

particle and the field particles are not equivalent since the field particles form a "bath". The field particles have a 
steady (given) distribution /(v') while the distribution of the test particle /(v,t) = NmP{v,t) changes with time. 
As a result, the Fokker-Planck equation derived by Hubbard [50| can be deduced from the kinetic equation ([55)1 
derived by Lenard [23 | and Balescu [25| and vice versa. 



F. The relaxation time of a test particle in a bath 

The derivation of the Fokkcr-Planck equation (j9ip . relying on a bath approximation, assumes that the distribution 
of the field particles is "frozen" so that /(v) does not change in time, at least on the timescale NtD corresponding 
to the relaxation time of the test particle (see below). This is always true for a thermal bath, corresponding to a 
distribution at statistical equilibrium (Boltzmann), because it does not evolve at all. For d > 1, due to collisions, 
any other distribution function relaxes towards the Boltzmann distribution on a timescale Nto. Therefore, for d > 1, 
only the Boltzmann distribution can form a bath (the other distributions change in time). The situation is different 
in d = 1 since a Vlasov stable distribution function does not change on a timescale of order Ntjj. Therefore, any 
Vlasov stable distribution function can be considered as "frozen" on this timescale and can form a bath. 
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Recalling that u ^ 1/N and f ^ N, the Fokker-Planck operator in Eqs. ([9T|) scales like 1/A''. Therefore, the 
distribution P(v,i) of the test particle relaxes towards the distribution /(v) of the bath [6^ on a typical time 

4"*'' r^Ntn. (93) 

This is the timescale controlling the relaxation of the test particle, i.e. the time needed by the test particle to 
acquire the distribution of the bath (for 3D stellar systems and plasmas, there are logarithmic corrections, so that 
tb^th _ (N/\nN)tD and tj?*'* ^ (A/lnA)i£, respectively). The tim escale tJi"*'' should not be confused with the 
timescale tR controlling the relaxation of the system as a whole (see Sec. Ill Cp . In d > 1 they have the same scaling but 
in d = 1 they are totally different since tji > Ntn (for spatially homogeneous systems). For one dimensional plasmas, 
the difference in behavior of "distinguished" particles that relaxes on a timescale At o and the overall population that 
relax on a timescale K^tn is shown in [s^, H^- As a clear sign of the inequivalence of these two descriptions (evolution 
of the system as a whole and relaxation of a test particle in a bath), we note that the Lenard-Balescu equation p6p 
conserves the energy while the Fokker-Planck equation (|9ip does not. 



IV. THERMAL BATH IN d DIMENSIONS 



A. Einstein relation 



We consider a thermal bath in d dimensions that is formed by field particles at statistical equilibrium having the 
isothermal (Maxwell-Boltzmann) distribution 



Substituting the identity 



f{-)=[^-^\'\e-y-\ (94) 



1^ = -/3m/(v)v, (95) 



in Eq. ([89|) . we obtain 



= -n{2nYm^p J dkdV k,j-^^^^^^S[k ■ (v - v')](k . v')/(v'). (96) 

Using the (5-function to replace k • v' by k • v, and comparing the resulting expression with Eq. (|87p. we find that 

Ff°' -A:j(v)/3toUj = -D\i{v)l3mvi, (97) 

where Z?|| (v) is the diffusion coefficient in the direction parallel to the velocity v of the test particles (see Eqs. (|109p 
and (|113p below that are valid for any isotropic distribution of the bath). The friction coefficient ^{v) = £'||(?;)/3m is 
given by an Einstein relation expressing the fluctuation-dissipation theorem. We stress that the Einstein relation is 
valid for the friction due to the polarization Fpoj, not for the total (dynamical) friction A that has a more complex 
expression due to the term djDij. We do not have this subtlety for the usual Brownian motion for which the diffusion 
coefficient is constant. For a thermal bath, using Eq. (|97|) . the Fokker-Planck equation ((53|) takes the form 



dP d 
dt dvi 



dP 



(98) 



where DjAy) is given by Eq. (|87p with Eq. (j94p . This equation is similar to the Kramers equation in Brownian 
theory [51| . However, in the present case, the diffusion coefficient is anisotropic and depends on the velocity of the 
test particle. For t — ^ -l-oo, the distribution of the test particle relaxes towards the Maxwell-Boltzmann distribution 
Pg(v) = (/3m/27r)''/2e"^""'/^ As we have seen, the relaxation time scales like t^"*'' Nto- 



B. General expression of the diffusion coefficient 



By introducing the representation 



^ e^-^-^— (99) 
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for the delta function in Eq. (IST]) . the diffusion tensor can be rewritten 



A,(v) = (27rf'^m 



dt / dhkikj 



u{kY 



|6(k,k.v)|^ 



(100) 



This equation can be directly obtained from the Kubo formula (|59|) . This shows that the auxiliary integration variable 
t in Eq. (jlOOp actually represents the time. Since the distribution function (|M)) is Gaussian, its three dimensional 
Fourier transform /(kt) is also Gaussian. The integration on time t, which is the one dimensional Fourier transform 
of a Gaussian, is therefore easily performed and we get 



A,(v)=.(2.)^('^ 



1/2 



pm 



kikj 



l{kf 



k |e(k,k-v)|2 



(k^ 



(101) 



Alternatively, this expression can be directly obtained from Eq. (|87p by introducing a cartesian system of coordinates 
for v' with the z axis taken along the direction of k and performing the integration. On the other hand, for the 
isothermal distribution (p4)) . the dielectric function is given by Eq. (jB4[) . Using Eq. (|B6p and noting that k- v is real, 
we can rewrite the diffusion tensor in the form 



A,(v) = / dkfc,fc,e-^^™(''-)^I?, ( ./^k. 



(102) 



where 



2(27r)" 



+ 00 



riikfk"^ 



[l-ri{k)B{x)Y +f^{kYC{xY 



■dk. 



(103) 



In the foregoing formulae, we have defined k = k/fc, B{x) = 1 — 2x6"^^ e^' dy, C{x) = ^xe^^' and rjik) = 
— {2'KYu{k)f3mp. In addition, «,„ = (l//3m)^/2 is the r.m.s. velocity in one direction and n = p/m the numerical 
density of particles. The function B{x) can be written B{x) = 1 — 2xD{x) where D{x) = e^^ dy is Dawson's 

integral. We note the asymptotic behaviors B{x) = 1 — 2x^ + ... for x — )- and B{x) ~ —l/{2x'^) for x — )■ ±cx). 
This function vanishes at x = 0.92414... and achieves a minimum at (1.502,-0.28475). Finally, we note that B{1) = 
-0.076159... and C(l) = 0.652049.... 

Equations (|102p - (|103p provide the general form of the diffusion tensor of a test particle in a thermal bath. The 



relaxation time can be estimated by v"^ 



hath 



^/D. Since the velocity of the test particles scales like {v ) ^ Dt, the 



relaxation time represents the time needed by the particle to acquire the typical velocity Vm of the bath. Alternatively, 
the relaxation time can be estimated by the inverse of the friction coefficient t^^*^ ~ 1/^- Using the Einstein relation 
^ ^ D(3m, we again obtain 4"*'' ^ Then, using D - vlJnR'^+^, which can be deduced from Eq. (jl03|) . and 

introducing the dynamical time t^*'' ^ R/vm, we obtain the scaling i^*'' ~ Ntu. A more precise discussion of the 
relaxation time will be given in Paper II. 

The foregoing expressions can be simplified for one-dimensional systems. In d = 1, the Fokker-Planck equation 
takes the form 



dP 
'dt 



d_ 

dv 



D{v) 



dP 

dv 



/3mPv 



(104) 



where the diffusion coefficient is given by 



Div) 



//3m 



(105) 



with 



Vi{x) = 



+00 



r]{k)^k 



2(27r)5 n [1 - r]{k)B{x)y + 7]{k)^C{x)^ 



■dk. 



(106) 
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C. Neglect of collective effects: Landau approximation 

If we neglect collective effects, which amounts to taking |e(k, k- v)| — 1 or, equivalently, B = C = in the equations 
of Sect. IIVBI we can write the diffusion tensor as 



(107) 



with 



2{2Try~i n Jo 



+ 00 



7]{kfk'^ dk 



(108) 



In the Landau approximation, the potential of interaction only appears in a multiplicative constant "Dd that fixes the 
relaxation time. Therefore, general expressions of the diffusion tensor can be given [s^ . 
In d = 3, introducing a spherical system of coordinates, we obtain 



D.,,{v) 



D\\{v)-]^D^{v) 



VjVi 1 
V 



(109) 



with 



G\ J^v 



erf ( ^/^v 



G{ 



(110) 



Giw) 



2 1 



t^e"* dt 



2w^ 



crilw) =e 



(111) 



where erf(w) is the error function. In the above expressions, £'||(f) and D±{v) are the diffusion coefficients in the 
directions parallel and perpendicular to the velocity of the test particle. We note the asymptotic behaviors 



In d = 2, introducing a polar system of coordinates, we obtain 



3, ^J^.{v) --+00 — 1=^ 

/3rn ^ 



(112) 



(113) 



with 



/3m 



(3m 



r , /^"^ 2 



where /„(a;) are the modified Bessel functions. We note the asymptotic behaviors 



T I ^"^ 



Z?||(0) =^2?2, D^{Q) = 'kV2, D\\{v)^+, 



2v^ 



Finally, in d = 1, 



(114) 



(115) 



D{v) = 22?ie-5 



(116) 
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D. Asymptotic behaviors and Debye-Hiickel approximation 

111 this subsection, we take collective effects into account but consider asymptotic expressions of the diffusion 
coefficient. 

For small velocities u — ^ 0, the function defined by Eq. (jl03p can be approximated by 

Therefore, when tj — ^ 0, the diffusion coefficient is given by Eq. (|107p with replaced by V^iQ). If we come back to 
the expression ()10ip of the diffusion coefficient, and use Eq. (|B4p . this amounts to replacing the dressed potential of 
interaction Udressed{k) u(fc)/|e(k, k • v)| by 

- l + (2.mfc)/?^p - 

When this formula is specialized to the Coulombian potential (27r)'^u(fc) = Sde^ /rn^k'^ , we see that Eq. (|118p is the 
Fourier transform of the Debye-Hiickel potential 



The replacement of the Coulombian potential by the Debye-Hiickel potential is a simple way to take into account 
collective effects in the kinetic theory by using the static results on screening. This procedure has been introduced 
phenomenologically by several authors in plasma physics before Lcnard and Balescu provided the rigorous solution 
(sec Paper II). The present considerations show that the Debye-Hiickel approximation is valid only for a thermal bath 
and for small velocities of the test particle. The results of Lcnard and Balescu arc more general because they account 
for dynamical screening and take the deformation of the polarization cloud, due to the motion of the test particle, 
into account. 

We now consider the case of large velocities v +oo. In d = 3 and d = 2 dimensions, it is not difficult to see that 
in the evaluation of D±(v), the function 'Ddix) can be replaced by 15^(0), while in the evaluation of D\\{v), it can be 
replaced by 

2(2^)-^-^ n Jo [l-j7(fc)S(l)]Vr;(fc)2C(l)2 

As a result, for v — >■ -|-oo, the diffusion coefficient D±{v) is given by Eq. (|107l) with Vd replaced by l?d(0), and the 
diffusion coefficient D^^{v) is given by Eq. (jl07p with Vd replaced by I?d(l). By contrast, in d = 1 dimension, there 
is no general expression for the large velocity behavior of the diffusion coefficient. It depends on the precise form of 
the potential of interaction. Analytical expressions will be given in Paper II for particular systems. 



V. OUT-OF-EQUILIBRIUM BATH IN d = 1 

In one dimension, we have seen that the Lenard-Balescu collision term vanishes due to the absence of resonance. 
Therefore, a stable steady state f{v) of the Vlasov equation does not evolve under the effect of collisions on a timcscale 
of order Nto- We can therefore develop a test particle approach in a bath characterized by an arbitrary Vlasov stable 
distribution /(u), not necessarily the MaxwcU-Boltzmann distribution of statistical equilibrium. Since the distribution 
f{v) does not change on the timescale t^*'' ~ Ntn corresponding to the relaxation of the test particle, we can consider 
that f{v) is "frozen" on this timescale, so the bath approach is self-consistent. 

In d = 1, the Fokker-Planck equation (|9ip takes the form 



dP 
'dt 



d_ 

dv 



p 



d\nf 

dv 



(121) 



with a diffusion coefficient 



D{v) — mf{v) 



\e{k,kv)\^ ■ 



(122) 
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According to Eq. (|B2|) of Appendix [B] one has 

\e{k,kv)\^ = 1 - 2nu{k)r J 



+00 ft 



fin) 



du 



+ ATT'^u{kf f [V) 



For the isothermal distribution (p4)) . we recover the resuhs of Section HVl For the water 
(IB3|) . we obtain 

D{v) = — I dk ^ ' 



2w„ 



1 



2-Ku{k)p 
2 — 2~ 



(123) 

distribution, using Eq. 

(124) 



for — w,„ < < '^n 
the form 



and D{v) ~ otherwise. When cohective effects are neglected, the diffusion coefficient reduces to 



D{v) — mfiv) 



+00 



4Tr^u{kfkdk. 



(125) 



It is proportional to the distribution function f{v). 

Equation (|121[) is similar to a Fokker-Planck equation describing the motion of a Brownian particle in a potential 
U{v) = — \nf{v) created by the field particles. The distribution function of the test particle P{v,t) relaxes towards 
the distribution of the bath f{v) for t — > +00. The typical timescale governing the approach of the test particle to the 
bath distribution is t^*'' ^ Nto- The fact that the distribution P{v,t) of a test particle relaxes towards any stable 
distribution f{v) of the field particles for c? = 1 explains why the distribution of the field particles does not change on 
a timescale ~ Nto- The situation is different in d > 1 32 1. 



VI. THE EFFECT OF AN EXTERNAL STOCHASTIC FORCING 

In a recent paper, Nardini et al. (57| have studied long-range interacting systems perturbed by external stochastic 
forces and derived a nonlinear Fokker-Planck equation describing the evolution of the system. We have independently 
considered the same problem. However, our approach is different and leads to a different final equation. The reason 
of this difference is unknown to us but we think that it is interesting to present our calculation in detail. We therefore 
consider the case of a system with long-range interactions that is submitted to an external stochastic potential $e(r, t)- 
The Klimontovich equation ^ remains valid provided that we make the substitution ^ + ^e- We shall treat 
the stochastic potential $e(r,t) as a small perturbation. We shall also assume that the system remains spatially 
homogeneous during the whole evolution. Then, decomposing fd = f + Sf and = $-|-i5<I' into a smooth component 
plus a fluctuating component, and making a quasilinear approximation as in Section [II B[ we find that Eqs. p2p and 
are replaced by 

^^^.^Sfm) + -^-{6fW^^), (126) 



dSf dSf df df „ 

at or ov ov 



(127) 



Equation (|126|) for the smooth component is exact while, in writing Eq. (jl27p for the fluctuation, we have neglected 
the nonlinear terms (5/(5$ — ((5/(5$) and (5/$e — ((5/$^). If we take the Fourier-Laplace transform of Eq. (|127p . we 
find that 



k • V — c 



(5$(k, w) + $e(k,w) 



'^/(k,v,0) 
i(k • V — w) 



(128) 



Substituting Eq. (|128p in Eq. (PH]) . we find that the Fourier-Laplace transform of the fluctuations of the potential is 
related to the initial condition and to the stochastic potential by 



(5$(k,a;) = (27r) 



d u{k) 
e(k,c.) 



dv 



(5/(k,v,0) 



i(k ■ V — uj) e(k, uj) 



1 - e(k,w) - 

^$e(k,Cj). 



(129) 



19 



We can use these expressions to compute the terms appearing on the right hand side of Eq. (|126|) . The first term 
can be written 

(SfVS^) = J dkj ^ J dk' J ^ik'e'('''-"*)e'('^'"-"'*)(,5/(k,v,tj)5l'(k',a;')). (130) 
Using Eq. (|128p . we obtain 

{Sf{k,v,uj)S^k',u')) = /'^ , \{S^k,u)6^k',u')) + (I.,(k,c.)5l.(k',c.')>l + (W'^y^'O)^'^'^('"^'^0) . (131) 



i{k • V — w) 



Some of the terms appearing in this expression have aheady been computed in Sect. IIIBI They lead to the Lenard- 
Balescu equation p6p . In this section, we shah only write down the new term arising from the stochastic forcing. 
Using Eq. (|129p . this term is found to be 



k • — 1 — f(k' ll)'] 

(./(k,v,.).$(k',.')) = j,) ^^-^'^'"^^-^'^'"^^^- ^'''^ 

We assume that the stochastic potential $e(r,<) is ^-correlated in time and that the spatial correlator only depends 
on the distance r — r'. Under these assumptions, we have 

(l'e(k, u;)de(k', u')) = P{k)6{k + k')6{uj + uj'), (133) 

where P{k) is the amplitude of the noise in Fourier space. Equation (IT32|) can then be rewritten 

k • ^ 1 — p(k I A* 

(<S/(k, V, u:)S^k', c.')) = , \J Pik)6ik + k')6{u; + lo'). (134) 

k • V — a; |e(k, aj)p 

The new contribution of the first term in the right hand side of Eq. (|126|) arising from the stochastic forcing is 
therefore 



(<5/V<5<i>), = -— I dk I — fc. ,_ J\ Pik)- (135) 



_ _ ^^fc^ ^-1^ l-6(k,C.) 

2t: J y 27r k • v — a; |e(k, a; 
Proceeding similarly, we find that the second term in Eq. (jl26p can be written 



du: ^ k-fi e{k,Loy 



Therefore, the net effect of the stochastic forcing is to lead to an anisotropic diffusion of the form 



Of 



/ forcing ^ i \ ^ 3 . 



with a diffusion tensor ^ jji ^ D^f where 

'■J ''J *j 



^^'t^l = -2^ J 2;^^-^-^ k.v-. HK^w ^^''^ 



Dm = -^fdkfp hk,^-^^^Pik). (139) 
■' 2tt J J 2tt k • V — w |e(k,a;)p 

The total diffusion tensor is 

Using the Landau prescription w — > w + iO^ and the Plemelj formula (PT|) . we can replace l/(k • v — w — iO+) by 
+i7r(5(k • V — io). Then, integrating over w, we obtain 

^,„....i^,^_L|,,^,^_m_. (141, 
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We note that, because of the presence of the dielectric function e(k, k • v) defined by Eq. ([22]) in the denominator 
of Eq. (|14ip . r'lj"'^'^*"^ [/] is a functional of f{v,t). For non-interacting particles, M(k) = 0, the diffusion coefficient 
arising from the forcing is given by Eq. (|14ip with e(k, k • v) ~ 1. This is just a constant depending on the amplitude 
P{k) of the noise. For interacting particles, the amplitude P{k) is "dressed" by the polarizaion cloud encapsulated in 
the dielectric function. The part Df^ of the diffusion tensor is explicitly computed in Appendix [U] 

Because of the stochastic force, the distribution function /(v,t) will diffuse and will not reach a steady state. 
However, a steady state can be reached if we assume that the particles, in addition to the stochastic force, also 
experience a friction force — due to the presence of an inert medium. If the friction force is small so that its 
infiucnce can be neglected in Eq. ()127|) . the calculation of the diffusion coefficient is not modified and we readily 
obtain a nonlinear Fokker-Planck equation of the form 



Of 



=^[D{r''V]^+tfv^]. (142) 

This Fokker-Planck equation is expected to relaxe, on a timescale of the order towards a non-Boltzmannian 
steady state determined, in a self-consistent manner, by the diffusion tensor _D^°^'^*"^ [/] . In fact, since this equation 
is nonlinear, it is not clear whether an equilibrium state will be finally achieved. 

Finally, if we develop a test particle approach and make a "bath approximation" , we can compute the diffusion 
tensor as in Section Hill In that case, we find that the new term in the Fokker-Planck equation (|?1]) due to the forcing 
and the dissipation is 

/ forcing+dissipation \ "^^i / 

where now /(v) is the time-independent distribution of the bath. In particular, for a thermal bath, using the expression 
(jB4p of the dielectric function, we get 

I3//''""^(v) = ^ / k^kjPik) __ ^^^^^ 



This Fokker-Planck equation relaxes, on a timescale of the order ^ ^. towards a non-Boltzmannian steady state 
determined by the diffusion tensor d(°^'^^^^ iy) that is itself determined by the distribution of the bath /(v). Simplified 
expressions will be given in Paper II. The results of this section can also be generalized to spatially inhomogeneous 
systems, using angle-action variables as in [ll| . 

Remark: We emphasize that the effect of an external forcing on the dynamics of systems with long-range interactions 
is physically different from considering Brownian particles with long-range interactions like in [s^,!!^,!!^- There exists, 
however, qualitative analogies. For Brownian particles, the effect of the thermal bath is represented, in the kinetic 
equation, by a linear Fokker-Planck operator of the Kramers type. The thermal bath breaks the conservation of energy 
and the system ultimately relaxes towards the Boltzmann distribution with the temperature of the bath, which is the 
canonical distribution. In the present case, the forcing and friction terms also break the conservation of energy (in 
a different manner) and the system ultimately relaxes towards a non-Boltzmannian distribution depending implicitly 
on the amplitude P{k) of the stochastic force in Fourier space. 



VII. CONCLUSION 



In this paper, we have reviewed and completed the existing literature on the kinetic theory of systems with long- 
range interactions. The Lenard-Balcscu equation introduced in the context of plasma physics has a wider scope, as 
it describes a large class of systems with (attractive or repulsive) weak long-range interactions. We have considered 
a test particle approach and given explicit expressions of the diffusion coefficient and friction force arising in the 
Fokker-Planck equation. We have stressed the analogies and the differences between the Lenard-Balescu equation 
describing the evolution of the system as a whole and the Fokker-Planck equation describing the relaxation of a 
test particle in a bath. We have also considered the effect of a stochastic forcing on the evolution of a system with 
long-range interactions and we have shown that it leads to an anisotropic diffusion with a diffusion tensor expressed 
as a functional of the distribution function itself. 

The general results that we have presented in this paper will be used, in Paper II, to investigate the kinetic theory 
of particular systems with long-range interactions. 
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Appendix A: Auto-correlation of the fluctuations of the one- particle density 

According to Eq. pT]) . we have 

/rlr Hr' 
' {Sf{r, V, 0)<5/(r', v', 0)> 

/rlr dr' , , 

(2^(2^^"'''''^' ' [(/4r,v,0)/4r',v,0)) - /(v)/(v')] . (Al) 

The expression of the discrete distribution function recalled in the first paragraph of Sec. Ill Bl leads to 
(/<i(r, V, 0)/<i(r', v', 0)) = ^ {5{y - r,)<5(v - v,)<5(r' - v,)5{w' - v,)> 

= m^Y^ (<5(r - r,)(5(v - v,)(5(r - r')(5(v - v')) 



+ m2 {S{r - r,)<5(v - v,)<5(r' - r,)S{V - v,)) 

= m/(v)J(r - r')<5(v - v') + /(v)/(v'), (A2) 

where we have assumed that, initially, there is no correlation between different particles (if there are correlations 
described by a smooth distribution /i2(r — r', v, v'), it can be shown that their contribution decays rapidly with time 
[23} so they have no effect on the final form of the collision term). Combining Eqs. (jAip and (|A2[) . we obtain 

((5/(k, V, 0)d7(k', v', 0)) = j^Sik + k')<5(v - v')m/(v). (A3) 

Appendix B: Dispersion relation 

We consider the linear dynamical stability of a spatially homogeneous steady state /(v) of the Vlasov equation ([8]). 
Decomposing the perturbation in a sum of plane waves of the form Sf ~ g«(k-r-(^t) ^^^^ substituting this decomposition 
in the linearized Vlasov equation (|13p we obtain the dispersion relation e(k,w) = 0, where the dielectric function is 
defined by Eq. (|22p . It determines the complex pulsation uj = Ur + it^i as a function of the wavenumber k. Taking 
the z-axis in the direction of k, we get 

e{k,co) = l-(2nfu{k) [ ^dv, (Bl) 



where we have noted v for Vz and f{v) for J f(v)dvxdvy. On the other hand, C is the Landau contour. Explicit 
expressions of the integral, depending on the sign of oJi, are given in [60| . In particular, when uii = 0, using the Plemelj 
formula ((3T|) . we get 



+ CXD r/ 



{k, u) = l- {27rf u{k)r I ' ^ dv - {27^fu{k)^7Tf' ('^) . (B2) 



Let us first consider particular distributions. 

• The waterbag distribution is defined by f{v) = p/2vm for —Vm < v < Vm and f{v) = otherwise (we note that 
v,n ~ corresponds to f{v) — pS{v)). Using f'{v) — {p/2vm)[S{v + Vm) — S{v — v^)], the dielectric function (|Bip is 
explicitly given by 

e(fc, w) = 1 + i27Tfu{k)p- ^7TT2 • (B3) 

The dispersion relation can be written = v'^k'^ + {2'K)'^u{k)k^ p. We note that a; is purely real or purely imaginary. 

• For the isothermal distribution (j94[) . introducing the notation ri{k) — -'(2Tr)'^u{k)/3mp, the dielectric function can 
be written 

eik, Lo) = l- r]{k)W (V^^) , (B4) 



22 



where 

is the M^-function of plasma physics [6lj. It has the exphcit expression 

2 2 2 

W{z)^l-ze^^j dy + i^-ze~^, (B6) 

for any z € C. We note that W^(0) = 1. On the other h and, fo r z e R, Wr(z) = B{z/V2) and ^^(z) C{z/V2) 
where the functions B{x) and C{x) are defined in Section HVBI 

Appendix C: The first component of the diffusion coefficient due to the forcing 

The first component (|138p of the diffusion coefficient due to the forcing can be written 

Using the Plemelj formula ([5T1) . we obtain 



^^^•'^ 47r 7 ' ^ |e(k,k-v)|2 ^ ' 47r2 7 7 ^ ^cj-k-v |e(k,cj)| 

On the other hand, the dielectric function (p2)) can be written 

e(k, w) = 1 + {2TTfu{k)r f ^' ^ dv - iTr{2Trf u{k) f k • ^,5(w - k • v) dv. (C3) 

J CO — K • V J OW 

Substituting Eq. JCH in Eq. wc get 

^^[^] - ^ /c^k.,fc,.(fc)7' jd.'^^^ (Rkj^ + RkjbyF) ^^'^ 
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